Dynamical generation of phase-squeezed states in a two-component Bose-Einstein 

condensates 
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As an "input" state of a linear (Mach-Zehnder or Ramsey) interferometer, phase-squeezed state 
proposed by Berry and Wiseman exhibits the best sensitivity approaching to the Heisenberg limit 
[Phys. Rev. Lett. 85, 5098 (2000)]. Similar with the Berry and Wiseman's state, we find that two 
kinds of phase-squeezed states can be generated dynamically with atomic Bose-Einstein condensates 
confined in a symmetric double-well potential, which shows the squeezing along spin operator S y 
and the anti-squeezing along S z , leading to sub-shot-noise phase estimation. 

PACS numbers: 03.75.Mn, 05.30.Jp, 42.50.Lc 



I. INTRODUCTION 

Recently, atom interferometer with a trapped Bose- 
Einstein condensates (BEC) has attracted much atten- 
tion due to its potential applications in quantum infor- 
mation and quantum metrology [H-0|- I n particular, a 
two-component BEC with the 'one-axis twisting' (OAT) 
interaction \S Z is widely believed to be a useful resource 
for preparing the OAT- type spin squeezed state @-[H|, 
as well as many-partite entangled state [12, [l3| . Using 
maximally entangled (i.e., the NOON) state rather than 
product coherent spin states (CSS), phase sensitivity of 
a linear (Mach-Zehnder, or equivalently Ramsey) inter- 
ferometer can be improved from the so-called shot-noise 
limit (SNL)A$ ~ 1/y/N to the Heisenberg limit (HL) 
~ 1/N [14H17| . where N is total particle number. How- 
ever, it is very difficult to create the NOON state. 

As a special case of spin squeezed state, Gaussian 
number-squeezed state (NSS) has been demonstrated 
in several experiments jl8l423j , which exhibits sub- 
Poissonian atom number distribution [3, [2l[ and long 
coherence time [22( due to the spin squeezing of relative- 
number operator S z . These features manifest it as a 
promisin g ca ndidate for the Heisenberg-limited metrol- 
ogy [24|, |25[. Dynamical generation of the NSS has 
been proposed in a two-component BEC by exploiting 
the interplay of nonlinear interaction and atomic tunnel- 
ing [26|, l27|. It was found that the variance of S z has 
an exact relationship with the visibility of the Ramsey 
signal (i.e., the phase coherence) [28l |29|. 

In this paper, we further investigate dynamical gener- 
ation of phase-squeezed states (PSS) with atomic BEC 
confined in a symmetric double well. We focus on a sym- 
metric input state of a linear interferometer, with which 
we derive easily the phase sensitivity and its relation with 
spin squeezing parameter. Two experimentally realizable 
schemes are proposed to generate the PSS. Similarly with 
the Berry and Wiseman's (BW) state [3(J, we find that 
it exhibits the squeezing of the spin operator S y and the 
anti-squeezing of S z . Numerically, we calculate power 
rules of the optimal squeezing with particle number N. 
Although the achievable squeezing is worse than that of 



the BW state, dynamically generated PSS are still useful 
to gain sub-shot-noise limited phase sensitivity [3l| . 

This paper is arranged as follows. In Sec. II we review 
quantum metrology protocol based upon the linear in- 
terferometer. For a symmetric input state, we derive the 
relation between phase sensitivity and the spin squeez- 
ing parameter [l4 |. In Sec. Ill, we investigate general 
features of the opti mal p hase-squeezed state proposed by 
Berry and Wiseman [30| ■ Two routes for creating the PSS 
are proposed with the two-component BECs, which uti- 
lize a cooperation between the nonlinear interaction and 
the Josephson tunneling of ultra-cold atoms confined in 
a symmetric double-well potential. Scaling rules of the 
optimal coupling and the resultant spin squeezing are 
obtained to compare with the BW state. Finally, we 
conclude in Sec. IV with main results of our work. 



II. A SYMMETRIC INPUT STATE OF LINEAR 
INTERFEROMETER 

Let us firstly review the achievable sensitivity of a lin- 
ear interferometer, which can be described by multiple 
rotations to the input state [32| . namely 
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where the three unitary transformations (from left to 
right) represent the output 50:50 beam splitter (or equiv- 
alently a 7r/2-pulse), the linear phase shifter, and the 
input beam splitter, respectively. In the Schwinger 
representation, the angular momentum operator S — 



(S x ,S y ,S z ) with 



jO L ), and S z 
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\{a\a R + a R a L ), S y = ^i(a' L a R 



L a,L — a R an), with denoting 



the annihilation operator for two Bosonic modes /i = L, 
R. Via a detection of S z (i.e., the light-intensity differ- 
ence or population imbalance) at the output ports, the 
dimensionless phase shift $ could be estimated with its 
precision quantified by 



(A5 Z ) C 



d(S z ) out /d$ 
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where (i) out and (Ai) out = [(i 2 ) out - (i) 2 ut ] 1/2 den ote 
the expectation value and the variance of an operator A 
with respect to the output state Bout- 
in terms of common eigenstates of S 2 and S z , \s,n) = 
(a[y+ n (4) s ~"|0)/V(s + n)!( S -n)! = \s + n, s - n) L , R , 
symmetric input state of a linear interferometer can be 
expanded as 
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where s — N/2 and the probability amplitudes c_„ = 
c n - It is easy to prove that the mean spin (S)i n — 
((S x )ini (Sy}ini (S z )in) is oriented in the x direction. Spe- 
cially, the zth spin component {S z )- m = due to c_„ = 
c n . In addition, expectation value of the ladder operator 
S 4 = (3 U - At . 



(£_)' = a' L an is given by (even N case) 



s-l 

E 



X_ n c„c* +1 = 2 X_„ Re [c„4 +1 ] 



n=0 



where X„ = (s + n) 1 / 2 (s — n + l) 1 / 2 , satisfying X Tn = 
X± n _|_i. Due to real (&+)&, we have (S y )i n = Ira(S+)i n 
i I and (S x )i n = Re(5_|_)i n ^ 0. Similarly, we can prove 
that (S x S y + S y S x )in = Im{S+)- m = and (S X S Z + 
Sj x ) in = Rc(S+(2S z + l)) in = 0. 

The above results, valid also for odd N case, en- 
able us to derive the phase sensitivity A$ in a sim- 
ple way. For instance, we get immediately the out- 
put signal (5' 2 )out = — (S x )i n sia.3> an d the variance 
(A^) ou t = [(AS z ) 2 n cos 2 $ + (AS^sin 2 *] 1 ^ J^, 
where we have used the relation: e l * SH (S' 2 ) fe e~ l * s « = 
(S z cos<I> — S x sin&) k with k = 1, 2. Inserting them into 
Eq. ([2]), we arrive at 



A$ = 



(A^) 2 n + (A^) 2 n tan 2 $ 



(S x ) 



which is minimized for the phase shift $ = 0, with the 
best sensitivity (A$) = (AS z )i n /\(S x )i n \. Here, the sub- 
script "0" denotes the achievable sensitivity for $ = 0. 
Following Wineland et al. [3], we introduce the squeez- 
ing parameter 



Cs 



(A$)p 
(A$) css 



N(AS Z 



\(S X ) 
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where (A<I>) CSS = 1/vTV, known as the shot-noise limit, 
represents the best sensitivity with product coherent spin 
states 331: 



of the mean spin \{S X }\ = s. Therefore, as an input of lin- 
ear interferometer, the coherent spin states lead to phase 
estimation limited by the shot noise 1/y/N (i.e., Cs = !)• 
How to improve the sensitivity beyond the SNL (i.e., 
Cs < 1) is one of the main subjects of quantum metrol- 
ogy [HI H3 • It has been found that several nonclassical 
states can reach the Heisenberg-limit (A$]n oc 1/N, such 
as the NOON state ^(|s,s) + \s,-s)) Qj-inj], and the 
twin-Fock state \s,0) (= \N/2, N/2) h>R for even N) [34| . 
However, experimental realization of these states is not 
an easy task. In this paper, we present experimentally 
available schemes for preparing the phase-squeezed state 
and its application in the quantum metrology. 



III. DYNAMICAL GENERATION OF 
PHASE-SQUEEZED STATE 

Previously, Berry and Wiseman have proposed an opti- 
mal phase-squeezed state |^)bw = S n c ™l s ' n ) with the 
probability amplitude [30j 



1 



2s + 2 



(6) 



Obviously, the BW state has a symmetric probability 
distribution about n = because of c_ n = c„, which in 
turn leads to the x-polarized mean spin (S) = ((S x ), 0, 0). 
As shown in Fig. [TJ one can also find that spin squeezing 
of the BW state is along j/-axis and the anti-squeezing 
along z-axis, i.e., 

(AS y ) 2 < s/2 and (AS^) 2 > s/2. (7) 

Here, the value s/2 (= N/4) denotes the standard quan- 
tum limit (SQL). 
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FIG. 1: (Color online) (a) Probability distribution jc n j 2 for 
the Berry and Wiseman's state, with c„ given by Eq. (|6]); (b) 
Quasi-probability distribution Q(9,(f>) — \{d, 4>\^}bw\ 2 on a 



Bloch sphere, for s — N/2 
\0, cj>) given by Eq. $9$. 



20 and the coherent spin state 



\s,±s) 
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Note that the mean spin of \s,±s) x is also parallel with 
the x-axis due to S x \s, ±s) x = In addi- 

tion, the two CSS satisfy minimal-uncertainty relation- 
ship [H: (AS Z ) 2 = (A5. y ) 2 = \{S x )\/2, with the length 



The phase-squeezed state is a useful resource to im- 
plement Heisenberg-limit metrology 30]. However, it 
should be pointed out that the PSS is not an input of 
linear interferometer, but a state after the first beam 
splitter. In other words, Eq. (UJ should be rewritten as 
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|*) out = e~ l % s *e l ' i ' s *\y} pss . Inserting it into Eq. ©, we 
obtain the best sensitivity (A$) = (AS , y ) pss /|(»5 , a ;)p SS | 
and hence the squeezing parameter £s = v^V(A$)o, 
where (S x } pss and (AS y ) pss denote the expectation value 
and the mean-square root of the variance with respect to 
|^) P ss- As an optimal PSS, the BW state gives rise to 
the Heinsenberg-limited sensitivity (A$)o ~ 2y/2/N and 
thus Cs ~ 2V2N- 1 I 2 for large enough particle number 

(N > 10 2 ) d|. 

In order to prepare the PSS, we now consider atomic 
Bose-Einstein condensate, which is tightly confined in a 
double well [2j, |22j, |36h38| | . By introducing two bosonic 
operators aum for left (right) mode-function </?l(r.) and 
the Schwinger's angular momentum operators as Eq. ([1}, 
quantum dynamics of the BEC system can be described 
by an effective spin Hamiltonian (with h = 1) [39rl44j ]: 



H 1 = 5S Z 



(8) 



where 5 denotes the potential bias, f2 the amplitude 
of Josephson coupling, and x the s-wave interaction 
strength. In the second term of the Hamiltonian © , an 
equatorial spin operator S 7 — S x cos 7 + S y sin 7 is intro- 
duced to simulate Josephson tunneling of the atoms with 
Rabi frequency fie 17 , where both the amplitude f2 and 
the phase 7 are tunable 0, 0] • The third term, arising 
from inter- and intraspecies nonlinear interaction, gives 
rise to the OAT-induced spin squeezing and many-partite 
entanglement 

Following original spin-squeezing scheme proposed by 
Kitagawa and Ueda [8(, and more recently by Riedel et 
al. |7j, we assume that a 7r/2-pulse is applied to the BEC 

to prepare the initial states Is, ±s) x = e~ l ^ Sy \s, ±s), 
which are coherent spin state [33j: 



\0, 



exp[—i9(S y cos0 — S x sin</>)]|s, s), 



(9) 



with polar angle 9 = tt/2 and azimuth angle cf> = or 
it. During the pulse, the two-mode functions </?l,r merge 
with each other such that S ~ x — an d hence the 



Hamiltonian H 



tt/2 



flSy, where the phase and the du- 



ration of the coupling are adjusted as 7 = flt w /2 = tt/2 
with Q » Nx [a, 0j- Next, the system evolves under 
the Hamiltonian Hq — QS X + X$ z f° r a duration t m - m . 
Nonzero x can be realized by splitting the two-mode BEC 
that allows for a suitable tunneling with the amplitude f2 
and the phase 7 = (or 7r), which is fulfilled by applying 
a second pulse Q. 

The interplay of the linear tunneling and the nonlin- 
ear interaction results in dynamical generation of Gaus- 
sian number-squeezed state I26M29I1 and Schrodinger's cat 
state (i.e., the NOON state) pUM50l |. Here, we show that 
it is also possible to create phase-squeezed state in the 
two-mode BEC. To understand how it works, we con- 
sider the Hamiltonian Hq with time-independent x an d 
il. For a fixed particle number N (— 2s), the system's 
energy is conserved, i.e., d(Ho)/dt — 0, which yields 



(S^(t)) — C — Sl{S x (t)) jx with a constant of integral 
C = s/2 ± fls/x for the initial CSS |s, ±s) x , that is 



(S 2 z (t)) 



1 ± 



2fl 

X. 



it 



{SS)) 



(10) 



where the upper sign corresponds to |s,+s)a; and the 
lower sign for \s,—s) x . It is interesting to note that the 
variance (AS" Z ) 2 = (§1) > s/2 only if the generated state 
is symmetric and evolved from Isj+s)^ (or |s,— s} x ) un- 
der Hq with > (< 0). Except the so-called 'over- 
squeezing' [51| . for which both (AS'z) 2 and (AS'y) 2 are 
larger than the SQL, the enhanced spin fluctuation of 
relative atom number implies the appearance of a PSS 
with (AS y ) 2 < s/2. For brevity, we will assume that 
both f2 and x are rea l an< i positive. 



A. Direct method to prepare the PSS 

Firstly, let us consider dynamical evolution of the ini- 
tial CSS Isj+s)^ under the system Hamiltonian Hq. To 
distinguish with previous scheme [47| , we adopt relatively 
larger coupling ratio f2/x ~ 1.1 JV (see blow Fig. 2]), which 
leads to a completely different dynamics due to the fact 
that the initial state, i.e., Eq. © with 8 = tt/2 and = 0, 
now corresponds to a stable fixed point in phase space. 
To confirm it, we plot time evolution of probability dis- 
tribution |c n (t)| 2 . As shown in Fig. H{a), the generated 
state shows symmetric probability distribution at 

any time t. For large N, the initial CSS \s, +s) x is a 
Gaussian namely 



M0)| 2 



1 

2? s 



2s 



1 



-n 2 /s 



(11) 



with its width determined by the variance AS Z = \J s/2 
(see blue line of Fig. HJb), also Ref. [H[). After a cer- 
tain duration i min ~ log 2 (iV)/(2iVx), the initial state 
will evolve into a phase-squeezed state, which can be in- 
ferred from a n in creased width of probability distribu- 
tion AS Z > y/s/2 (red line of Fig. [2]) . Such a result can 
be explained from Eq. (|10|). Indeed, we can prove that 
the mean spin of \*f>(t)} aligns with the x-axis and its 
length (S x (t)) < s. Therefore, from Eq. (jTU]) we have 
(A5* 2 ) 2 > s/2 and hence (AS y ) 2 < s/2, except for the 
over-squeezing. 

The appearance of the PSS can be illuminated 
by calculating quasi-probability distribution Q{6 1 <j)) = 
I (6», <j)\y(t)) | 2 on the Bloch sphere 0, where the CSS \6, <t>) 
is defined by Eq. (JSJ , with the polar angle < 9 < tt and 
the azimuth angle < <fi < 2ir. As depicted in the left 
plot of Fig. [3£a), an elliptic shape of the Q function is 
found for the spin state \^{t)) at t min , which shows the 
squeezing along S y and the anti-squeezing along S z , sim- 
ilar with that of the BW state. In Fig.[3Jb), we plot time 
evolution of the squeezing parameter (solid red line) 



4 



and the fidelity between the prepared state with the BW 
state, F(t) = |bw(*I*(*))| 2 (dotted line). For instance 
at t — 0, the fidelity is given by 



F(0) ~ \ exp 



1 



1 



(12) 



where we used Eq. (j6|) and the approximate result of 
Eq. (JTTJ). For the case N = 2s = 40 and fi/x = 1.134AT, 
our analytical result predicts F(0) ~ 0.675, coincident 
with exact numerical simulation 0.672. From Fig. [3jb), 
we find that the generated PSS at time t min — 6.95 x 
10 _2 x _1 shows a high fidelity F = 0.982 and a local 
minimum of the squeezing parameter Cs = 0.507 (—2.949 
dB), which is the optimal squeezing that this scheme can 
reach for N = 40. 
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FIG. 2: (Color online) (a) Probability distribution |c ?l | 2 of 
the spin state evolved from s, +s) x as a function of time t (in 
units of x^ 1 ); (b) Snapshot of \c n \ 2 for the initial state (blue 
line) and the PSS at time i m i n = 6.95 x 10 _2 x _1 (red line with 
points). Other parameters: fi/x — 1.134JV and N — 2s — 40. 



It is necessary to investigate the optimal coupling ra- 
tio fi/x and the minimal value of Cs (i- e -i the maximal 
squeezing) for any finite N. As an example, we con- 
sider the exactly solvable case with N = 2. Analytical 
results can be obtained as the following: (S x ) = 1 — 
2fc*\ „„j t a cg2 = i[i_^ s i n 2(a)], where £ = 
(x/2) 2 , denoting half of level-spacing between the 



^sin^a) and (AS y ) 



ground state and the second-excited state. Therefore, we 
get the maximal squeezing Cg = (Ax 2 + 1)/ (2x + 1) 2 with 
x = fi/Xi which occurs at times tf. = (k + ^)tt/£ for an 
integer k = 0, 1, • • ■ . Minimizing £| with respect to x, 
we further obtain the optimal coupling ratio fi/x = 1/2 
and thus the best squeezing Cs = 1/V2- For three- 
particle case, (S x ) = |[1 — fjsm 2 (£t)] and (AS y ) 2 = 



3 p_ _ 2 X (n- X ) sin 2( ft )] ) where £ = ^2 _ n)( + x 2 

maximal squeezing also occurs at times with 



cl 



(x 2 -x+ \){x 2 -3x + 3) 
x 2 (x — l) 2 



The 



(13) 



from which we find that the optimal squeezing Cs = 0.762 
is attainable for the coupling ratio x = fi/x = 2.794 
(~ 0.931 AT with N = 3). Counter-intuitively, a slightly 
weaker spin squeezing is obtained in a comparison with 



previous N = 2 case. This is because the generated PSS 
in the two-particle system is also maximally-entangled 
NOON state, which results in the Heisenberg-limited sen- 
sitivity (A$)o = 1/N and the optimal squeezing Cs = 
1/VN (with N = 2). 

The optimal squeezing for A^ > 3 has to be determined 
numerically. Due to computation limit, only finite par- 
ticle number 2 < N < 200 are shown in Fig. [4] Our 
numerical simulations show that the optimal squeezing 
with Cs ~ A^ - 21 can be obtained for the coupling ratio 
fi/x ~ 1.1 AT (see balls of Fig. [4j. Although the achiev- 
able squeezing is worse than that of the BW state (empty 
circles) [35| , it is still useful to achieve the sub-shot noise 
in the phase estimation due to Cs < 1- 
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FIG. 3: (Color online) (a) Quasi-probability distribution 
Q(8, (j>) of the dynamical generated PSS (left), the NSS (mid- 
dle), and the PSS: e i(,r/2) ^ |*} nss (right); (b) The squeezing 
parameter Cs (solid red) and the Fidelity F(t) (dotted) as 
a function of time t (in units of x _1 ) f° r the initial states 
\s, +s} x with s = N/2 = 20 and fi/x = 1-134AT. The dashed 
blue line corresponds to evolution of Cs from the initial CSS 
s, -s) x , with s = 20 and fi/x = 2.664. 



B. Indirect method to prepare the PSS via a 
rotation of the NSS 



Next, we consider dynamical evolution of another ini- 
tial CSS \s, —s) x under governed by the Hamiltonian Hq. 
According to Law et al. 26[, it is possible to generate 



dynamically a Gaussian number-squeezed state, which 
shows the squeezing along S z and the anti-squeezing 
along S y (see middle plot of Fig. [2Ja)). During time 
evolution, the generated state \^(t)) is always symmetric 
around n — 0. More specially, the probability amplitudes 
obey C- n (t) — c n (t) for even A" case, or c_„(i) = — c n (t) 
for odd N case, leading to the re-polarized mean spin [271 ]. 
Previously, the optimal coupling ratio fi/x as a func- 
tion of A" has been obtained by minimizing the reduced 
variance over the SQL: = 2(AS Z ) 2 /s, which was re- 
ferred as the number-squeezing parameter 0, [2^, |46| 
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and is measurable by a detection of the phase coherence 
= \(S x )\/s. For the symmetric state, Eq. (fTU|_gives 
us an exact relationship between and [2& l29|. 

To confirm whether the prepared NSS is useful for 
metrology or not, we have to investigate the squeezing 
parameter £s = £n/<7^ [20]. As shown by dashed blue 
line of Fig. ^[b), one can find that the prepared state 
\^(t)) at time t m j n = 0.159x _1 shows the optimal squeez- 
ing Cs = 0.322 (about -5dB) for the case N = 40 and 
0,/x — 2.664. In Fig. 21 we investigate power rules of 
the optimal coupling ratio fi/x and the best squeezing 
Cs- Our numerical results (the red crosses) show that 
fl/x ~ 0.917V 27 and Cs ~ 1.24Y- - 37 . Such a scaling is 
worse than that of the BW state Cs ~ A -0 - 5 . However, 
as an input state of linear interferometer, it is sufficient 
to achieve the sub-shot-noise limited metrology. 
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FIG. 4: (Color online) (a) The optimal coupling ratio fi/x as 
a function of N for the PSS (balls) and the NSS (red crosses), 
which fit with 1.1 JV (thin solid) and 0.91147V 2704 (thick blue 
line), respectively; (b) Power rule of the optimal squeezing Cs 
for the BW state (empty circles), the PSS (balls), and the 
NSS (red crosses). For large N (> 10 2 ), Cs of the BW state 
and the NSS fit well with 2^/2N' 1/2 (thin solid line) and 
1.2432Y -0 - 3715 (thick blue line), respectively. 



Before closing, it should be mentioned that the Gaus- 
sian numb er- squeezed state has been prepared in several 
experiments |18h23| , which exhibits sub-Poissonian atom 
number distribution [lH [2l[ and relatively longer coher- 



ence time compared with that of the product CSS [22j . 
These features manifest it as a promising candidate for 
quantum metrology [13, [25[ . Recently, Uys and Meystre 
have found that the NSS gives rise to the Heisenberg- 
limited phase sensitivity [2J| (see also Ref. [25|]). If a 
Gaussian NSS is fed into the interferometer as Eq. (JTJ, 

the action of the input beam splitter e l ^/ 2 ' Sx , equiva- 
lent with a 7r/2-pulse, rotates the prepared NSS about 
negative s-axis by an angle 7r/2, which results in a 

PSS: |#} pss = e^^^l^nss, as shown by right plot of 
Fig. Eta). This is another method to generate phase- 
squeezed state [3|, |4| . 

IV. CONCLUSION 

In summary, we have demonstrated that it is possible 
to generate phase-squeezed state in a two-mode Bose- 
Einstein condensate. Considering a symmetric input 
state of linear interferometer, we show that it is easy 
to derive the best sensitivity and its relation with the 
squeezing parameter Cs- As an example, we have dis- 
cussed general features of the Berry and Wiseman's state: 
the x-polarized mean spin and the Heisenberg-limited 
phase sensitivity. Two experimentally realizable schemes 
are proposed to generate the phase-squeezed states. We 
find that it can be prepared from an initial coherent 
spin state \s 1 +s) x without any rotation. Alternatively, 
starting with \s, — s) x , a number-squeezed state is pre- 
pared and then transformed into phase-squeezed state 
via a 7r/2-pulse. The optimal coupling and the maximal 
squeezing are determined based upon a wide range of nu- 
merical simulations. Our results show that the generated 
phase-squeezed states can achieve sub-shot-noise-limited 
quantum metrology. 
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